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Abstract 

It is argued that a diffusion may be ergodic even though the drift field has un- 
bounded outward-directed parts. The discussion employs stochastic and numerical 
methods. 
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1 Introduction 

We shall investigate the ergodicity of the following stochastic differential 
equation (SDE) 

' dX t = b(X t )dt + dW t , 

Xo = (, U 

where ( 6 M 2 and (Wt)t>o is a 2-dimensional Brownian motion. The drift 
b : M 2 — » M 2 is given by 

b x = -Ax\ + 12xiX2, b 2 = -I2xlx 2 + 4^2- (2) 
Here a stochastic process X = (X t )t>o is termed ergodic if there exists a 
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Figure 1: The drift field b according to (TSJ) 



unique invariant probability measure /i for X and if for every / 6 L 1 (R d ,/x) 
it holds that 



where denotes the law of X with starting point (. Condition ([3]) is also 
called the law of large numbers of the process X. Under very general assump- 
tions it can be shown that the existence of an invariant probability measure 
is equivalent to X being positive recurrent, see 0, Theorem 4.1]. As we shall 
see it is difficult to classify the diffusion X, which solves equation (pQ), as be- 
ing transient, recurrent or positive recurrent. This is due to the funny drift 
b itself. Before proceeding we fix some notation. 

Notation: Ef[-] shall denote the expectation with respect to the law Pf. 
The standard inner product of two vectors x,x' G M. d is written as x ■ x'. 
We set r = r(x) := y/x ■ x. Furthermore, B(M. d ) is the space of bounded, 
real- valued measurable functions and C 2 (M. d ) stands for the space of twice 
continuously differentiable functions. The indicator function of a set A will 
be denoted by 1a- 

Let us have a closer look at b defined in ([2]). Let e r := {xx,X2)/r be the 




P c -a.s., for all C e R d , 



(3) 



radial unit vector field, then we may distinguish two cases 

(a) e r ■ b < 0, if \x 2 \ < \xi\, 

(b) e r ■ b > 0, if \x 2 \ > \xi\. 

Therefore, in case (a) we have an inward drift and in case (b) an outward 
drift. Moreover, evaluating e r ■ b with x\ = rcosip, x 2 = rsimp, <p G [— 7r,7r], 
we see that in both cases there results a growth behavior proportional to 
r 3 and that the proportionality constant is positive and increases as \(p\ — > 
n / 2 for {7r/4 < \(p\ < 37r/4} and decreases, being negative, as \cp\ — > 
or \<f\ — > Ti for {\(p\ < 7r/4} and {\(p\ > 37r/4}, respectively It equals for 
\ip\ = 7r/2, \(p\ = 3tt/A. This means that moving away from the origin there are 
areas with unbounded inward drift as well with unbounded outward 

drift. Moreover, it should be noted that curl(&i,6i) = — 12x\x 2 which shows 
that it is impossible to write the drift as b = — VT^ for some differentiate 
function V. 

The paper is organized as follows. In section[2]we collect some background 
material on the different large time asymptotics that a stochastic process 
may show. These properties are discussed in terms of the coefficients of the 
generator of X. In section [3] we return to the investigation of the diffusion 
solving ([1]). It is made evident that the latter diffusion cannot be handled 
using the standard criteria of section [2j Finally, the numerical results for the 
long-term behavior of ([3]) are presented which give strong indications that 
the diffusion solving (jTJ) is ergodic. 



2 Criteria for transience and recurrence 

Here we discuss criteria that allow for a complete classification of a diffusion, 
at least if the latter are applicable. To this end we consider a class of more 
general SDEs, namely 

dX t = b(X,)dt + dW„ 

x = C, (5) 

where (Wt)t>o is a d- dimensional Brownian motion and b : IR d — > IR d is a 
smooth vector field. The semigroup (S t )t>o, S t : B(M d ) — > B(M. d ), associated 
with X is given by S t f(C) := [/(^) !{<<»?}] > where r\ denotes the explosion 
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time of X. We shall assume that the generator L of X uniquely determines 
(S t ) t > and that on C 2 (R d ) C D(L) it is given by 

L = |A + 6V, with formal adjoint L*(-) = div((|V - &)(•)). 

Since a priori we do not demand the semigroup to be strongly continuous the 
generator L might be given also in the weak sense, see e.g. [U Section 2.2] 
and (21 Section 2.3]. From now on processes sharing the foregoing properties 
will simply be called diffusions. 

In [1] Bhattacharya has given quite general sufficient criteria in terms 
of the coefficients bi that allow to decide whether a diffusion is transient, 
recurrent or that it even admits a finite invariant measure. We recall that 
the diffusion X is said to be recurrent if for any ( 6 l rf and any open set 
U C M. d we have 

P^(X tn G U for a sequence of t n 's increasing to oo) = 1. 

If for any ( G R d it holds that 

P ? (limr(X;) = oo) = l, 

t— >oo 

then X is called a transient process. Finally, X is said to be positive recur- 
rent, if E^[<7{/] < oo, for all ( G M. d and all open sets U C W 1 , where o\j is the 
first entrance time of U. 

Adopting the same notation as in [1] , we introduce the following functions 

d 

C(x) = 2 x jbj(x) 
i=i 

P{r) = sup (d - 1 + C(x)), (3{r) = inf (d - 1 + C{x)) (6) 

\x\=r ~ \ x \= r 

7(r) = rm iu , 1{r)= rm du 

Jro u j ro U 

With these preparations we may quote the following two theorems from 

PQ 

Theorem 2.1 (a) If for some r > 

exp{— I{u)}du = oo, (7) 
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then the diffusion with generator L is recurrent, 
(b) If for some r > 

POO 

/ exp{— I_(u)}du < oo, (8) 

Jro 

then the diffusion with generator L is transient. 

Theorem 2.2 (a) The diffusion with generator L is recurrent and admits a 
finite invariant measure (unique up to a constant multiple) if there exists a 
r > such that (T7p holds and 

poo 

/ exp{I(u)}du < oo. (9) 

J rn 



' ro 

(b) If there exists a ro > such that |?[) holds and 

-N / rs \ / pN 

N->co 



lim / exp{— I(s)\ ( / exp{I(u)}du J ds / / exp{—I_(u)}du = oo, 

t-^Jro ' \Jr J I Jr 

(10) 

then the recurrent diffusion does not admit a finite invariant measure. 

Henceforth a unique finite invariant measure will always be normalized 
to become a probability measure. In the following examples we explore some 
possible scenarios. For simplicity only radially symmetric problems are con- 
sidered and we set J(r) := J(r) = I(r). 

Example 2.1 (a) Let b = in then X t = W t and an elementary cal- 
culation gives by and (TJ]) that Brownian motion is recurrent in d = 1,2 
and transient for d > 3. 

(b) Let V{x) = —r~ a and b = — VV, mt/i a > 0. In t/iis case we find 

I (d-l-2au a ) du = -lnr ) + 2(r~ Q -rn"). 

elementary calculations one deduces that (0) /io/ds /or d = 1,2, implying 
recurrence, and that ra/zd /or d > 3 entailing transience. On the other 

hand, evaluating the quotient of the l.h.s in / TiOj) gives the following estimate 

quotient of integrals of ffTUl) > C(r )N, 



5 



for N large enough, with C{tq) some positive constant. This shows that there 
is no finite invariant measure for any dimension. Still there is a a-finite 
invariant measure v{dx) = e~ 2V dx for every dimension, since the stationary 
Fokker-Planck equation has a solution 

L*(e- 2V ) = div((§V - b)e- 2V ) = 0. 

(c) Let V(x) = r a and b = — W, where a > 0. Here we obtain 

I(r) = jf ( ^ d ~^ ~ 2au a -^ du = In/" 1 - 2r a - lnr^ 1 + 2r°, 

which entails properties (f?p and Now the diffusion is recurrent with finite 
invariant measure. 

Remark. The functions V in examples (b) and (c) are not smooth at the 
origin. But the latter may be easily rendered such. 

It is clear that the diffusive part dW t tends to spread a "particle" moving 
according to (JSJ). This spreading may be compensated by a sufficiently strong 
drift b. As can be seen from the previous example (c), if the driving field b 
is inward directed and the gradient of a potential that grows like a positive 
power, then we are guaranteed an invariant probability measure \i on M. d . In 
fact in the case of example (c) the diffusion is even ergodic, meaning that for 
every / G L 1 (M d , jj) we have 

lim - / f(X t )dt= I fdfji, P c - a.s, for all ( G R d . (11) 

T^oo T J Q J R d 

Notice that, for / bounded and measurable, property (fTTi) immediately im- 
plies 

lim i / StfiQdt = [ fdfi, C 6 M d . (12) 

T^oo 1 J J Rd 

A large class of processes for which properties (fTTT) and ( 1121) hold is given 
by recurrent diffusions with invariant probability measure that satisfy in 
addition the following conditions (cf. [H Section 1.3]) 

Condition (A) The semigroup (S t )t>o maps Co(M d ) into 60(1^) and is 
strongly continuous. 

Condition (B) There is a Borel measure v on M. d with support on M. d 
and a strictly positive function p t (x,y) continuous for (t,x,y) G (0, 00) x 
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IRL x M. d such that the transition probability can be written as K t (x,dy) 
p t (x,y)v(dy). 



3 An ergodic diffusion with strong outward 
drift 

We now come to the discussion of the diffusion (X t ) t > which solves ([T]). 
The peculiar form of ([2]) is related to a complex differentiation, namely if 
z = x\ + i%2 and g — (x% + ix^f then 

bl = _ Re f*( g) \ h = _i m 

\dz J \dz 

We know at least that X has continuous sample paths, maps B(M, d ) into 
C(M ) and is conservative, i.e. its explosion time is almost surely infinite. 
The first property is a consequence of our assumptions, see e.g. [5]. The 
second and third properties are consequences of Lemma 2.5 and condition 
(4.1) of P. 

Let us turn to the transience and recurrence properties of our diffusion. 
The absolute maxima and minima of (1 + C(x)) given the constraint \x\ = r 
are found to be /3(r) = 1 + 8r 4 and (3(r) = 1 — 8r 4 , respectively. This gives 
I(r) = lnr — lnro + 2(r 4 — r 4 ) and J(r) = lnr — lnro — 2(r 4 — r 4 ), hence 

/ exp{-7(u)}du = r e 2r o / -e'^du < oo (13) 

J ro J ro H 

poo poo -t 

/ exp{-I(u)}du = r e~ 2r o / -e 2u4 du = oo. (14) 

J ro J ro " 

This shows that none of the criteria ([7]) and applies. As far as the 
existence of a finite invariant measure is concerned we use the estimate 

f N I f N 1 4 

quotient of integrals of (JTUD < C(r ) / (s — r )ds / / -e 2u4 du, (15) 

Jr / Jro u 

for some positive constant C(r ). Since the r.h.s of f|T5l) tends to as N — > oo 
criterion ( flOj) is not valid. Hence, even if the diffusion were recurrent we could 
not tell whether it possesses an invariant probability measure. 
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In order to get some insight into the problem of the diffusion ([T]) several 
numerical tests have been performed whose aim was to study the behavior 
of 



1 ' T 



h ■= f J f(X t )dt, (16) 

as T increases. As fx is an expression depending on the sample paths a 
strong Taylor approximation (of order 1.5) was used for the latter. More 
precisely, the strong Taylor scheme is implemented through (cf. [7J p. 355, 
eq. 4.12]) 

Yn+i = Yn + + AW k + \L% k A 2 + L%AZ\ k = 1, 2, (17) 

where 

AW k = U k VA, AZ k = \A^ 2 (U k + ^U k ), 

with Ui , C/| being independent standard Gaussian random variables and 
where A is the time step. Moreover, 

L° = b 1 1 + b 2 d 2 + \ [d\ + dj) , L 1 = dx, L 2 = d 2 . 

The corresponding results are shown in Figure [2] with A = 10~ 4 . Note that 
for / = 1a, A a Borel set, equality (J3J) entails 

hm 1 / l A (X t )dt = f x(A). (18) 

Hence indicator functions comprise a natural class of test functions to be 
used. A glance at Figure [2] clearly displays a convergent behavior of expres- 
sion ffTB]) for different starting points and different positions of the indicator 
functions. One also sees that the limits decay rather fast and that the dif- 
fusion takes ever more time to stabilize the farther away we move from the 
origin. On the other hand, if convergence of (Tl8|) is observed then the limit \i 
is a natural candidate for an invariant probability measure. Note that accord- 
ing to [21 Theorem 2.2.5] all transition probabilities K t (x, dy) of the diffusion 
X are equivalent which by Doob's theorem, see [3, Theorem 4.2.1], gives 
the uniqueness of the invariant measure /i. This means that the numerical 
findings support the assertion on ergodicity of X. 
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Figure 2: Behavior of /t, with / = l Bl (x — x ) where B\ = {x E M 2 \ r(x) < 
1} and for x the following values hav^been chosen, (a) (0, 0), (b) (2, 0), (c) 
(0,2), (d) (3,0), (e) and (f) (0,3), (gj 2(1/^,1/^), (h) 3(1/^,1/^). 
The starting points ( have been picked randomly from [—10, 10] 2 . 
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